Modelling Washboard Road : 
from experimental measurements to linear stability analysis 
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When submitted to the repeated passages of vehicles unpaved roads made of sand or gravel can 
develop a ripply pattern known as washboard or corrugated road. We propose a stability analysis 
based on experimental measurements of the force acting on a blade (or plow) dragged on a circular 
sand track and show that a linear model is sufficient to describe the instability near onset. The 
relation between the trajectory of the plow and the profile of the sand bed left after its passage is 
studied experimentally. The various terms in the expression of the lift force created by the flow of 
granular material on the plow are determined up to first order by imposing a sinusoidal trajectory to 
the blade on an initially flat sand bed, as well as by imposing a horizontal trajectory on an initially 
rippled sand bed. Our model recovers all the previously observed features of washboard road and 
accurately predicts the most unstable wavelength near onset as well as the critical velocity for the 
instability. 

PACS numbers: 45.70.-n, 81.05.Rjel, 62.40.+i 



I. INTRODUCTION 

When submitted to the repeated passages of vehicles, 
gravel and sand roads can develop a ripply pattern known 
as washboard or corrugated road. Not only is the bum- 
piness of the track a disturbance to drivers, but it also 
causes a loss of adherence and control and is therefore a 
real hazard. A similar phenomenon occurs on train, tram- 
way or metro tracks and is known as rail corrugation. 
Due to its obvious economic importance, rail corruga- 
tion, which is caused by wear or plastic deformation, has 
been extensively studied as reported in the engineering 
literature [l[ through experiments , field work Q and 
theoretical and numerical analysis The formation 

of a washboard road over a sand bed was first studied 
experimentally in |lOI|. and theoretical models have been 
proposed in [IJ, [l2[ but include an ad-hoc diffusion term 
whose physical meaning remains unclear. Recently, it has 
been shown that a washboard road can develop when a 
simple inclined blade or plow (instead of a rolling wheel) 
is dragged along a sand road [H, [l4| . The formation of 
an instability during a unique passage of a plow under 
its own weight over an initially flat surface has also been 
studied in viscous or viscoplastic fluids [HI, [l6[ . 

The aim of the present paper is to derive a model for 
the washboard instability caused by the repeated pas- 
sages of a plow over a sand bed, from experimental mea- 



surements of the forces acting on the plow. Previous work 
has focused on force measurements in the simpler case of 
a blade at constant altitude steadily dragging a mound 
of sand over a flat sand bed [17[ • Herein, we extend these 
results by using force measurements to probe the mecha- 
nical response to an oscillatory excitation. These measu- 
rements are then used as the basis for a linear stability 
analysis of the washboard road instability. 

The paper is organised as follows : Section|TT]gives a list 
of previous experimental observations and results that a 
model should recover. Section llTIl presents the experimen- 
tal methods while section IIVI introduces the framework 
and assumptions of the linear stability analysis. The rela- 
tion between the trajectory of the plow and the shape left 
in the sand bed after its passage is discussed in section IVl 
The expression of the lift forces acting on the plow is dis- 
cussed in section fVTl Finally the predictions of the linear 
stability analysis are compared to experimental results in 
section IVIIl 



II. PREVIOUS RESULTS 

Our previous work has established a number of charac- 
teristic features which a model for washboard road ought 
to reproduce, (i) There exists a critical velocity, v c , be- 
low which the sand bed remains flat when perturbed and 
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above which any irregularity will develop into a regular 
rippled pattern, (ii) This critical velocity increases with 
increasing mass of the plow or wheel, following a power 
law. (iii) For velocities greater than the critical velocity 
the rippled washboard road pattern appears only gra- 
dually over several passages of the plow. While the initial 
wavelength, A keeps a finite and well-defined value, the 
amplitude of the ripples grows continuously from zero to 
a saturated value. 

In the present paper we focus on the onset of the in- 
stability, i.e. on velocities greater yet close to the critical 
velocity. The saturation and coarsening of the pattern 
remain to be studied in further detail. 

Previous work has shown that in the case of wash- 
board road caused by a plow, compaction of the sand 
bed has not been observed nor seems to play any role 14| , 
contrary to the case of a rotating wheel. This may be due 
to the fact that when using a plow the granular material 
is constantly shuffled at every passage of the plow. There- 
fore compaction will be neglected throughout the present 
stability analysis. Note that this assumption would be 
highly debatable had a wheel been used. 

For the range of velocities and masses of the plow used 
we found that the " horizontal" length of the plowed ma- 
terial, Lo on figure (typically a few centimeters), re- 
mains smaller than the wavelength of the washboard road 
pattern, A typically ten centimeters. Its influence is brie- 
fly discussed in section IVI CI but unless otherwise men- 
tioned it will be neglected. 




Figure 1. Sketch of the experimental setup. A plow is dragged 
at a constant horizontal velocity on a 5-m-long circular track 
filled with sand. 



III. EXPERIMENTAL SETUP 

The experimental setup consists in a circular track of 
length 5 m (average diameter of 80 cm), width 25 cm, 
filled with a 20-cm high layer of granular material. The 
granular material used is sand-blasting sand of size ran- 
ging from 200 fim to 400 /im. The plow consisting of a 
flat blade (inclined at 45°) of width 14.5 cm is dragged at 
a constant horizontal velocity v, ranging from 0.1 ms _1 
to 2 ms _1 . These values of the speed are high enough 
to produce a continuous-flow regime and low enough to 
avoid a gaseous regime. The plow is either attached to an 
arm whose end is free to rotate (figure [5^) or mounted 
on a translation stage (figure [5b)- 

In the former case (figured), the length of the arm is 
50 cm while the amplitude of the ripples near onset is ty- 
pically of the order of a few millimeter. Therefore, we will 
consider that the arm remains nearly horizontal and the 
pendulum motion of the arm is neglected, so that any 
change in either the horizontal speed or in the inclina- 
tion of the plow due to a change in altitude is neglected. 
The vertical dynamics of the plow is therefore governed 
only by its own weight and the lift force caused by the 
plowed material, while its horizontal motion is imposed 
by a motor and the washboard road pattern may deve- 
lop. Compared to most engineering work [l[ our system 
is simplified in that the plow has no tire nor suspen- 
sion. Movies in the supplementary material show both 
the steady-state washboard instability on real-time and 
the growth of the ripples through stroboscopic images. 
The vertical position of the plow is recorded using a ma- 
gnetic angle sensor (ASM-PRAS1) placed on one end of 
the arm. Unless otherwise mentioned, prior to any ex- 
periment the sand bed is made flat by dragging (over 
several tens of rotations) a gradually rising vertical blade 
around the track. 

In the second case (figure [2b) the plow is rigidly moun- 
ted on a vertical translation stage (5 /im accuracy), while 
its horizontal velocity remains imposed. Two force sen- 
sors (Testwell KD40S) are used to measure the vertical 
lift force acting on the plow. In this second case the tra- 
jectory of the plow is imposed through a NIUSB6259 
card. The sand bed profiles prior to and after the pas- 
sage of the plow are measured using two laser telemeters 
(optoNCDT 1302 from Micro- Epsilon not shown on fi- 
gure [2Jd) of 0.02 mm accuracy. Obviously no washboard 
road pattern can appear using this setup used only to 
probe the mechanical response of the material to an im- 
posed trajectory. 

The following variables are defined : x is the horizontal 
position of the lower tip of the plow (set to zero at time 
t — 0), y(x — v(t — (n — l)? 1 )) its vertical position where 
n is the number of passages and T the duration of one 
passage, h n (x) is the profile of the sand bed after the n th 
passage of the plow (ho being therefore the initial profile 
prior to the experiment) (see figure [5]). 

Note that since the horizontal velocity is imposed there 
always exists a correspondence between the position x 
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and the time t. Any variable can therefore be differentia- 
ted with respect to x or t, although for convenience h n 
will be expressed as a function of x and y as a function 
oft. 

The sand bed being largely thicker than the amplitude 
of the ripples near onset it can be considered infinitely 
deep (having halved the thickness of the sand bed in our 
experiments showed no noticeable effects). The reference 
of the vertical variable y and h n is therefore arbitrary. 
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Figure 2. Sketch of the experimental setups, a) Plow atta- 
ched to an articulated arm and free to move vertically, b) Plow 
mounted on force sensors attached to a computer-controlled 
translation stage which allows one to impose the trajectory. 



IV. FRAMEWORK OF THE LINEAR 
STABILITY ANALYSIS 

The washboard road pattern emerges from the interac- 
tion between the plow and the sand bed. The trajectory 
of the plow is governed by Newton's second law where 
the vertical lift force, fi, acting on the plow may depend 
on the plowed mass M, the vertical position of the plow 
y, and the sand bed profile h and their derivatives. Simi- 
larly the shape of the sand bed depends on the trajectory 
of the plow, while the new profile left after the n th pas- 
sage, h n (x), may depend on the lift, the previous profile 
h n -i(x) and the trajectory of the plow y[x). 

The plow position y is governed by its dynamics 
(weight and lift force) while an erosion/deposition law 
reshapes the bed profile h (compaction being neglected). 
Although these two processes are coupled and simulta- 
neous a convenient way to represent the coupling is the 
schematics of figure [3] commonly used in the engineering 
literature [H H, Q . In the following section the reshaping 
of the sand bed is discussed whereas the dynamics of the 
plow is studied in section IVII 
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Figure 3. Diagram representing the coupling between the 
dynamics of the plow and the reshaping of the sand bed due 
to erosion and deposition. 



As mentioned in the introduction the present paper 
proposes a linear stability analysis near the onset of the 
washboard road instability. It will be shown that al- 
though nonlinearities may be needed to account for the 
saturation and coarsening of the ripple pattern, a linear 
model accurately reproduces the experimental data. The 
linear response is probed using sinusoidal excitations (of 
either the vertical position of the plow, or of the initial 
bed thickness) of wave-number k, or corresponding an- 
gular velocity w = vk. As mentioned above, the origin of 
the vertical axes is arbitrary and does not play any role. 
Since compaction is here neglected the average height 
and thickness profile is constant and can be set to zero. 
Therefore, the system can be described using complex 
variables M(t) = (M) +M e lwt , where (M) is the plo- 
wed mass on time-average, and similarly y = A y e tLJt and 
h = A h e lkx . As a reminder, based on previous experi- 
mental observations the four following assumptions are 
made, (i) A 2-dimensional model is suitable, (ii) com- 
paction is neglected and (iii) the length of the plowed 
material (Lq on figure [5]) is smaller than the wavelength 
of the pattern and (iv) near the onset of the instability 
studied here the plow remains in contact with the sand 
bed. 



V. RESHAPING OF THE SAND BED 

In rail corrugation, the time evolution of the bed profile 
can involve wear or plastic deformation [l[ and is gover- 
ned by the normal and tangential forces acting between 
the wheel and the track, as well as by the trajectory of the 
wheel y(t) and the previous shape of the track, h n ^i(x), 
these three variables being coupled. Here in the case of a 
plow running over a sand bed the sand is simply eroded 
and redeposited (compaction being neglected). 

In order to investigate the relation between h n (x) and 
the other variables the plow is dragged at constant ho- 
rizontal velocity v while a vertical sinusoidal trajectory 
of angular velocity, w, (and corresponding wavenumber 
k = uj/v) is imposed to the plow mounted on the trans- 
lation stage : y(t) — A v e iut . The profile, h n (x) is then 
recorded. Figure |4] shows a typical example of the out- 
come for v = 0.5 ms _1 , A y — 0.6 mm and k = 21 m -1 
(corresponding wavelength A = 30 cm). 

The first important result is that there is no phase-lag 
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simply given by : 
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Figure 4. a) Plots of the trajectory of the plow y (solid 
symbols) and the profile left after the passage of the plow, 
h n (open symbols), for k = 21 m^ 1 and v = 0.5 ms -1 . The 
profile h n measured after the passage of the plow is translated 
by the distance to the tip of the plow, b) and c) Ratio between 
the amplitude of the profile h n and the trajectory, \ = h n /y 
as a function of k (for v = 0.5 ms -1 ) and v (for k = 21 m _1 ). 



between the imposed vertical position, y, and the profile 
left after the passage of the plow, h n (see Figure 2k). 
Our experiments have shown that this remains true for 
the range of values of v, A y and k of interest in this paper 
(near the onset of the instability, v is ranging from 0.2 to 
2.0 ms -1 and k is ranging from 10 to 125 m" 1 ). 

Another major result is that the amplitude of h n is 
slightly smaller than that of y, meaning that the plow 
leaves a footprint in the sand which is similar yet not 
identical to its trajectory. The ratio between the two am- 
plitudes, x, is plotted as a function of wavenumber and 
velocity (in figure 0]b and c). The ratio x is clearly in- 
dependent of the wavenumber and depends very weakly 
on the velocity. The ratio x can be seen as an imprinting 
efficiency ratio. It seems intuitive that it should tend to- 
ward one for low velocities since when trying to carve a 
shape in a sand bed, one should be as delicate (and slow) 
as possible. Note however that over the range of velocities 
of interest x remains rather constant at x ~ 0.9 and any 
dependence on v will be neglected in the following. Mo- 
reover, the ratio x is also found to be independent of the 
amplitude of the imposed trajectory, A v , for the range of 
small amplitudes studied (A y < 5 mm). In the present 
case where the trajectory is imposed by the translation 
stage the profile left after the n th passage of the plow, h n , 
is found to be independent of that prior to this passage, 
h n —\. In the following, we will therefore consider that 
the profile left after the n th passage of the plow, h n (x) is 



K(x) = x y(x) 



(1) 



where x = 0.9±0.05. Although this value is close to unity 
it will emerge that the difference between the trajectory 
of the plow and the shape left in the sand bed is a key 
element of the washboard road instability. 



VI. DYNAMICS OF THE PLOW 
A. Jerk equation 

In a situation where the plow is free to move vertically 
its dynamics is simply governed by the lift force acting 
on the plow // (its horizontal position being imposed) : 



my 



-mg + fi 



(2) 



Our recent study [17[ has shown that in a uniform mo- 
tion over a flat bed (constant horizontal velocity and no 
vertical motion), the lift force is simply proportional to 
the plowed mass M and is independent of the velocity : 
fi = fiMg, where p, is an effective solid friction coefficient 
which depends solely on the inclination of the plow (or 
angle of attack), and p — 0.56 for an inclination of 45°. 
The goal of this section is to extend our previous analysis 
to non-uniform motion over a rippled sand bed. 

From mass conservation (and having neglected the 
length Lq, see figure [2]) the rate at which mass is gai- 
ned or lost is simply given by : 



M = pv (h n -i(x) - h n (x)) 



(3) 

bulk density 



where p is the 2D density of the sand (p - 
of sand x packing fraction x plow width) . 

In a steady state (horizontal or sinusoidal motion) 
equations and f3| simply give a relation between the 
average plowed mass and the mass of the plow : 



fi(M) = m 



(4) 



On an initially flat sand bed (h n _i(x) — 0) differentia- 
ting equation ([2]) and using equations ([T]) and ([3]) yields : 



h n (x) + K 3 h n (x) = 



(5) 



where K — (/j.gpx/imv 2 )) 1 ^ . This equation known as 
a jerk equation governs the dynamics of third-order os- 
cillators and charged particules in motion in their own 
electromagnetic field. Equation ([5]) has oscillatory solu- 
tions (of wavenumber y/3K/2) but with an exponentially 
growing amplitude, which is in total contradiction with 
the experimental observation of a critical velocity, at least 
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in the range of masses and velocities under study. There- 
fore additional dissipative terms are needed in equation 
([5|). Note however that the wavelength predicted by the 
jerk equation §5§ is rather close to that observed near 
the onset of the instability (i.e. the most unstable mode 
for a velocity close to the critical velocity). Indeed, for 
v = 0.8 ms _1 and m = 0.25 kg, equation ([5]) predicts a 
wavelength A = j= = 0.41 m whereas the experimental 



wavelength of the pattern is \ exp — 0.35 m. 

To the first order described in this section, the lift 
force depends on the plowed mass M, i.e. on J h n and 
j h n -\ as seen from equation [3l The simplest mathe- 
matical form is to include the derivatives of h n and 
h n -i, whose physical meaning will be discussed below : 

/; = /(/fc n ,//in-i,An)ftn-ijftn>^n-i)- Under the li- 
nearity assumption used here it will be assumed that 
the effects of all further terms are additive. Section fVI Bl 
is devoted to the study of the case h n -i = whereas 
section IVI CI will focus on the case where h n — and 
h n _ 1 {x) = A h e ikx . 



B. Lift force over an initially flat bed 

This section aims at determining the dependence of 
the lift force on h n and h' n in the case of an initially flat 
and horizontal bed, h n -\ = 0. The plow is mounted on 
the translation stage and the vertical force is recorded as 
the vertical position is imposed. 

An initial mass of sand (M ) is plowed by lowering the 
plow in the sand bed. From this position (at t = 0) a 
sinusoidal motion is imposed to the plow : y — A y e lujt = 
A y e lkx . The sand bed is flattened at every rotation using 
a vertical blade attached to another arm which erases any 
pre-existing profile. This allows the data to be averaged 
over typically 30 oscillation periods. The lift force reads : 



In order to measure the coefficients b and c and to 
study their dependence on the average mass (M) and 
velocity v, the transfert function H_ is defined : H_ = 
f { /y- Equation ^ yields the following band-stop filter 
expression : 
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(7) 
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Figure 5. a) Absolute value and argument of the transfer 
function \H\ for v = 0.3 ms" 1 and (M) = 0.58 kg. The solid 
lines correspond to equation (0) the dashed line corresponds 
to tj c for which H_ is real and \H_\ is minimal, b) lo c as a 
function of v showing a linear dependence for (M) — 0.58 kg. 



f l (t)=fiMg-by-cy (6) 

where b and c are real and positive coefficients whose de- 
pendence on the average plowed mass (M) and on the 
velocity is key to the stability analysis. Note that in the 
steady regime studied here equation © could be expres- 
sed as a function of h n and h' n using equation ((T|). 

Although the two additional terms by and cy are in- 
troduced as the simplest linear extension of the lift force, 
their physical meaning is easily understood. The first ad- 
ditional term can be seen as a restoring force acting on 
an intruder, adding an extra force which is proportional 
to the penetration depth y. Clearly when pushing an in- 
truder into the sand bed (y < 0), the additional force is 
positive, hence the negative sign in —by. Such linear res- 
toring forces have been reported in the literature [18l - l21| . 
The second term may be seen as a contribution of the ver- 
tical penetration speed (or equivalently of the local slope 
of the profile when ^ 0). Again the additional force 
should be positive when penetrating into the sand bed 
(y < 0), hence the negative sign in — cy. 



Equation ([7]) predicts that there exists a minimum in 
\ H \ when H_ is real and negative (ip — arg(g) — ir), i.e. 

1/2 

for u! = Ld c = {^9PX V I C ) ■ Therefore the coefficient c 
can be computed from the value of w c while the coef- 
ficient b is simply the minimum reached by \H\. In the 
framework of equation ((6|) measuring min(|_ff|) and uj c is 
sufficient to experimentally determine the values of the 
two coefficients b and c. The experimental data of the 
absolute value has a higher noise level than the argu- 
ment and hence, the value of w c was measured from the 
argument ip{ui c ) = ir (dashed line on figure [5]). 

The transfer function H_ was measured for various plo- 
wed masses, (M), and velocities, v , and a typical example 
is shown on Figure [S] (v — 0.3 ms -1 and (A/) = 0.58 
kg). The experimental data shows the features predic- 
ted by equation ([7]) (minimum in \H\, reached for H_ real 
((/> = 7r), ip ranging from tt/2 to 37r/2) and both the ab- 
solute value and the argument are well fitted by equation 
([7]) (solid lines). For the small amplitudes imposed here 
(typically A y < 5 mm) we found no dependence on the 
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excitation amplitude A y . These results validate the linear 
form of equation ©. 

The values of lo c computed from the transfer function 
are plotted in figure [SJd as a function of v for an ave- 
raged plowed mass (M) — 0.58 kg . It appears that 
lo c is a linear function of the velocity v : uj c oc v + v , 
where vq = 0.1 ± 0.01 ms -1 . The physical meaning of 
this dependence remains to be explained. Still this em- 
pirical expression along with measurements of min(|_ff|) 
for various plowed masses and velocities allows to plot 
the two coefficients b and c as functions of (M) and v 
(figure [5]) • The following observations emerge : over the 
range of parameters of interest, b does not show any sys- 
tematic dependence on (M) and is a linear function of 
v (Fig. [B^ and |B]d) and the following expression will be 
used b = Bq v, with Bq = 560 kgm _1 .s _1 . The coefficient 
c appears to be proportional to the mass (M) (solid line 
of slope one in figure [5]:) whereas its dependence on v 
has already been determined from the empirical expres- 
sion of u) c (solid line on figure [5]) . Overall the following 
expression emerges : c = Cqijl(M) g p\v / (v + vq) 2 , where 
Co is a constant Cq = 8.7 10 -2 m.kg -1 . 



C. Lift force over a rippled sand bed 

The previous section has established empirical expres- 
sions for the contribution of the derivatives of y(t) (or 
equivalently h n (x)) to the lift force, In this section 
the contribution to // of a pre-existing profile prior to the 
n th passage of the plow is studied by shaping a sinusoidal 
profile (h n _i — Ahe lkx — Ahe lut ) over which a horizon- 
tal trajectory is imposed (y = 0, and hence h n = 0) with 
a blade plowing an average mass of sand (M) . 

Similar additional terms to the lift forces are expected 
although with an opposite sign. Indeed a bump in the 
sand bed (ft.„-i > 0) will create an extra positive force 
as will a positive slope h' n _ 1 : 



400 



fl = fiMg + b h n -i + c h n - 



(8) 



where b and c and coefficients whose dependence on the 
averaged plowed mass (M) and on the velocity v has to 
be determined. 

It is expected that the role of y in equation ([6]) and 
/i„_i in equation (J5J) are symmetrical and we propose 
that b — b and c — c. A new transfer function can be 
defined as if' = / /h n _ 1 and is plotted as a function 
of u) on figure [7] (for (Al) — 1 kg and v — 1.5 ms -1 , 
one can note that v is high compared to the previous 
experiment, this is to prevent the effect of the length Lq 
and make sure that A ^> Lq). The solid lines correspond 
to the predictions of equation (|8]) using the expressions 
of b and c obtained in section IVI Bl They show a good 
agreement with the experimental data and validate the 
proposed expression of the lift force. 
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Figure 6. (a),(c), Coefficients of the additional terms b and 
c as a function of (M) (for v = 0.5 ms -1 ) and v (for (M) = 
1 kg), b appears to be proportional to v (solid line of slope 
one) and independent of M. c is proportional to M and its 
dependence on v is given by the expression of uu c found in 
figure (solid line). 
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Figure 7. Absolute value and argument of the transfer func- 
tion for M — 1 kg and v = 1.5 ms -1 . The solid lines 
correspond to equation ([8]) with parameters inferred from 
section IVI Bl showing good agreement with the experimen- 
tal data. 



VII. STABILITY ANALYSIS 

A. Exponential growth rate 

Having empirically established the mass and velocity 
dependence of additional terms in the lift force in sec- 
tion I VII using a plow whose trajectory is imposed, the 
equation of motion of a plow free to move vertically can 
be given by combining equations ([!]), © and ([S]) under 
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the assumption of additive forces made above : 

my = -mg+(j,gM(t)+b(v) (h n -i-y)+c((M),v)(h n -i-y) 

(9) 

Differentiating equation ([UJ and using equation (fTJ) the 
growth rate of the instablitiy, a = In | ^ | , can be ex- 
pressed : 

1 / (a - 7 fc 2 ) 2 + {/3kf \ 
a = ln 7 — " , ' 9 , Q , , qN9 + In x (10 

where a = /ig p / '(mn 2 ) , /3 = B a /(mv) and 7 = 

Co#Px/( u + u o) 2 - 

Experimentally this growth rate is measured by pre- 
paring an initially rippled pattern (similarly to sec- 
tion IVI CI) and instead of imposing a horizontal trajec- 
tory to the plow it is left free to move vertically. After 
several rotations of the plow the ripples either vasnish or 
increase depending on the wavelength of the pattern and 
on the velocity of the plow. Figure [8] shows two typical 
examples of these behaviors. During the first rotations 
(at least 10), the amplitude of the ripples can be fitted 
to an exponential curve and the growth rate is calculated 
from the exponental fit. 
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Figure 8. h n as a function of the number of rotations for 
v = 1.1 ms" 1 and A = 195 mm (□) and v = 0.8 ms -1 and 
A = 216 mm (o). Straight lines are exponential fits of the data 
over the first ten rotations. 

Using the values of the coefficients B and Co determi- 
ned experimentally in the previous sections, the growth 
rate a can be computed from equation [TU] as a function of 
the wavenumber k (or corresponding wavelength A) for 
various values of the average mass (M) and velocity v. 



Figure [9J (top) shows plots of a for (M) = 0.3 kg and 
for velocities ranging from 0.4 ms -1 to 1.2 ms -1 . For 
comparison, the experimental growth rate per rotationis 
shwon in figure [UJ bottom. 
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Figure 9. (Top) Predictions of equation , growth rate a 
as a function of A for velocities ranging from 0.4 ms -1 to 
1.2 ms~ from bottom to top : (110 j) . (Bottom) Experimental 
measurements for v = 0.6 ms -1 (x), v = 0.8 ms -1 (•) and 
v — 1.1 ms -1 (o). Above a critical velocity v c there exist po- 
sitive values of a corresponding to the onset of the instability. 

The predictions of the model are in good qualitative 
and quantitative agreement with experimental results 
and show an important feature : for low velocities, the 
growth rate remains negative for all wavelengths, mea- 
ning that any perturbation will be gradually eroded, whe- 
reas for high velocities there exists a positive maximum 
corresponding to the fastest growing mode or the most 
unstable wavelength. 

B. Critical velocity and fastest growing mode 

The first velocity for which there exists a positive value 
of the growth rate, a, is the critical velocity of the instabi- 
lity, v c . The model derived not only recovers the existence 
of a critical velocity but also predicts the fastest growing 
mode : the theoretical and experimental most unstable 
wavelengths are in excellent agreement as shown in fi- 
gure [51 

It should be noted that there is some degree of un- 
certainty as to the asymptotic behavior of a since the 
expressions for the additional terms in equation ([6]) were 
determined over a limited range of wavelengths and ve- 
locities near the transition. 

The critical velocity v c for which the first positive va- 
lue of a is observed depends only on one parameter : 
the average plowed mass (M), i.e. the mass of the plow 
m = fi(l\t). Figure [TU] (top) shows the values of v c in- 
ferred from our model for a plow mass ranging from 100 
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Figure 10. (Top) Critical velocity v c as a function of the mass 
of the plow m (o), the straight line corresponds to the predic- 
tion of the model and follows a power law with an exponent 
0.35. (Bottom) Experimental wavenumber as a function of the 



predicted one. The straight line shows k e 



g to 2 kg. The critical velocity follows a power-law as a 



function of this mass 



,0.35 



, as found in previous 



experiments [i"4j . yet with a lower exponent (v c oc m 25 ) 
Although it is satisfactory to recover a power-law, the ori- 
gin of the discrepancy in the exponent remains unclear. 
However the relevance of a power-law over only one de- 
cade is debatable. Again expanding the predictions of 
equation (| 10[) to a wider range of masses, m, would be 
questionable since the expression of the lift force was de- 
termined in a somewhat narrow range of parameters near 
the onset of the instability. 

Finally the value of the predicted most unstable wa- 
venumber, kth, can be compared to that measured expe- 
rimentally, k exp . The latter was not measured from the 
plots of a but instead was simply deduced from the pro- 
file of the sand bed after a few rotations. These values of 
kth and k exp depend on the mass of the plow m as well as 
on the velocity, v (clearly for v > v c (m)). For simplicity, 
they are plotted not as a function of m and v but instead 



k exp is plotted as a function of kth on figure [TOl (bottom) . 
The agreement between the predicted and experimentally 
measured most unstable wavenumber is excellent for all 
values of the mass and velocity as shown by the solid line 
of slope one. 

VIII. CONCLUSION 



We have presented a linear stability analysis based on 
experimental measurements of the lift force acting on a 
blade plowing a mound of sand on a sand bed. From pre- 
vious experimental work |17j | which focused on the case 
of a steady plow over a fiat surface an equation for the 
dynamics of the plow was derived. However we found 
that further terms were needed to recover the observed 
features of the instability. The empirical expressions of 
these additional terms were obtained by probing the me- 
chanical response of the system to sinusoidal excitations, 
both in the trajectory of the plow and in the initial profile 
of the sand bed. From this the growth rate of a sinusoidal 
perturbation was calculated and we showed that above a 
critical velocity, which depends on the mass of the plow, it 
displays a positive maximum corresponding to the most 
unstable or fastest growing mode. The critical velocity 
as well as the most unstable wavenumber were compu- 
ted and showed excellent agreement with experimental 
measurements . 

The expression of the lift force was obtained for a range 
of parameters near the onset of the washboard road in- 
stability. Hence the model recovers all the experimental 
observations near the threshold of the instability but fails 
to include the saturation of the pattern observed in expe- 
riments. Further measurements would be necessary and a 
thorough modelling will require to include nonlinearities. 

Finally it would be interesting to extend the present 
study to the case of a rolling wheel. However this remains 
challenging since on the one hand the transport mecha- 
nism is clearly different and may involve plastic defor- 
mation and on the other hand compaction is expected to 
play a crucial role. 
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